


«s£» 

mmk 

IStew^^ 

^Soi'y^S 


£>?'fXfir 


I ... .' 


Lfia£rii> 


^NATIONAllftEGHNI€^M 

^NR0RMATIG^iSER^|E^ 

Sr Department, bliCommarea 


$£#1 








lil 



REMTECH INCORPORATED 


RTR 008-4 


FOREWORD 

This report presents results of work performed at 
REMTECH' s Huntsville office. The work was conducted for 

Marshall Space Flight Center (MSFC) In response to require- 

/ ' ’ 

merits of Contract NAS8-28115. 
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Section 1 
INTRODUCTION 

Several technical areas were encompassed in providing support for Booster 
thermal environment test work. These areas Included cavity flow heating, rare- 
fied flow heating and Impulse operated model research and testing. Cavity 
flow heating problems were studied with respect to the proposed altitude 
control motors for the Space Shuttle. Available literature on this subject 
was reviewed and analytical predictive methods were surrmarlzed for use in 
planning testing work. Rarefied flow heating data was reviewed and correlated. 
The study showed the importance of considering rarefied flow conditions In 
launch thermal environment prediction. Impulse operated model research 
and testing was conducted to provide a basis for understanding and designing 
such models for booster thermal environment testing. During the course of 
this work technical reports were completed which document the various tech- 
nical faclts encompassed. This set of reports. Refs. 1 to 13, provides 
comprehensive documentation and thus a synoptic presentation method is used 
herein. The most salient features of the work are given in the following 
sections. 
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Section 2 
CAVITY HEATING 

Propulsive attitude control systems which are intermittently used 
during hypersonic atmospheric flight are discussed and analyzed as an open 
cavity flow heat transfer problem in Ref. 1. The literature which presents 
open cavity flow characteristics is reviewed' and discussed In terms of ACS 
nozzle-cavity heating. A description of existing analytical models for 

i ...... 

prediction of heat transfer within cavities Is given In Ref. 1 and are 
summarized here. All of the models presented were developed for two-dimen- 
sional laminar flow (see Fig. 1). Thus, unless three-dimensional effects 
within ACS nozzle cavities are found to be negligible these models will 
require modification to be applicable to ACS nozzle flow. Furthermore, 
other models for turbulent flow will be required for anticipated flight 
conditions. ? 

2.1 BURGGRAF'S MODEL > 7 ; - 

Burggraf (Ref. 14) modeled the steady separated flow in two-dimensional 
rectangular cavities at high Reynolds numbers. The cavity flow was considered 
to consist of a single constant enthalpy, rotating Inviscid eddy which produces 
a thin boundary layer along the cavity walls. The laminar momentum and energy 
equations were linearized and the Prandtl number assumed to be one to obtain 
a tractable solution. The heat transfer distribution was normalized to using 
Chapman's average value transferred across the dividing streamline. The 
expression relating the preseparation heating to the distribution within 
the cavity is given in Table I. 


2 



RTR 008-4 


REMTECH INCORPORATED 

Burggraf also examined the limitations of his analysis with respect to 
the assumption of an In viscid central core. The model is no longer valid 
when the boundary layer thickness becomes equal to one half of the cavity 
width, W, for cavities with H>W, This condition is related to a critical 
Reynolds number given by Eq. (2) of Table 1. Thus for the theory to be 
applicable the local edge Reynolds number, Eq. (1) of Table 1 must be 
greater than the critical Reynolds number given by Eq. (2) of Table 1 . 

j As part of the present investigation, Burggraf' s model was programmed 
to provide a predictive method. Due to the simplicity of the heat transfer 
correlation, the program was used to parametrically evaluate the heat trans- 
fer distribution. The results of this parametric analysis are given in Fig. 

2 for W/H values spanning the range of expected ACS nozzle cavities geome- 
tries. . v i: \ ' r .• 

2.2 HODGSON'S MODEL 

Hodgson (Ref. 15) visualized cavity flow as determined by an oncoming 
boundary which separates and forms a shear layer which in turn produces an 
equivalent wedge flow down the reattachment wall. The equations corres- 
ponding to this conceptual model are given in Table 2 and are discussed 
below. 

The effect of finite upstream boundary layer was Introduced through 
its influence on the dividing streamline velocity, U d , as determined from 
results of Denison and Baum (Ref. 16). In the present work, the theoretical 
curve was least squares fit for numerical computation purposes. 

Cavity fluid properties are computed using an average cavity tempera- 
ture, Eq* (4), and the boundary layer edge pressure. The dividing streamline 
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temperature, T d , Is computed using a modified Crocco's relationship, [Eq. (3)]. 

One of the significant contributions of this model Is the compressibility 
correction given for the reattachment zone length, L, expressed in Eqs. (6) 
and (7)* This length Is used in an expression for the local ve! col ty along 
the reattachemnt wall which was developed for Incompressible flow, [Eq. (8)]. 
Reattachment wall heat transfer is then computed using heat transfer relations 
[Eqs. (9) to (13)] for a flat plate having a varying freestream velocity. 

| The theories of Burggraf and Hodgson provides the most easily used methods 
for the engineer to obtain an estimate of the heating distribution within a 
cavity. Several additional methods and aspects of this problem are presented 
and discussed in Ref. 1. 

The cavity heating information given in Ref. 1 was applied to a specific 
problem in Ref. 2. This reference documents a study of the thermal environ- 
ments for two proposed cavity regions of the External Tank. Thermal environ- 
ments for the exit region cf the LOX Tank Vent Line and exit region of the 
Interbank Deorbit Retro were computed for two trajectory points corresponding 
to peak heating trajectory points. Heating rates for specific points within 
the cavities and near the cavities on the external surface are given In tabular 
form in Ref. 2. ' 


4 
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Table 1 

Burggraf: Laminar 2-D Cavity Heating Equations 


1. Re, - PeUeX ° 
° 

Nomenclature 

U e - Boundary layer (BL) edge velocity 
p e - BL edge density 

y e - BL edge viscosity 

/Y \ / \ 

X Q - BL running length to cavity 

2. Re = 240(— ) M + H] 

y - Distance along cavity reattachment 

cr \W/ V w / 

wall from external surface 


^ps“ ^separation convective heating rate 
q(y) -Local* heating rate 

3. If Re v > Re then 
x 0 cr 

W,H- See Fig. 1 



4. Rleman zeta function (a rational approximation) 
n-0 

a = + 0.803323 

9 

* x s * 3.89728 
a. = + 2.55002 

i ... 

a * - 1.19121 

. S - 

a « + 0.308284 ■ 

k 

a f - 0.0335024 

5 i 
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Table 2 

Hodgson: Laminar 2-0 Cavity Heating Equations 


1, $* * W/X Q (Flat Plate) 


U, 


2. ip = f(S*) from Fig. 5 of Denison and Baum (Ref. 16) 
u e 


Least Squares Fit 


Uj 

— * a + a s + a c 2 + a £ 9 + a c* + a 


3. 


where € = 


log (S*) 
10 . 


a * + .51706 

0 

a = + .088707 

1 

a = - .049822 

2 

a * - .0020116 

3 

a = + .0047757 
«* 

a = + .0006549 
s 


W 


d l/T 


aW T w) 


U. 2 


4. T c = .5(T d + T M ) 


5. Re u = 

" w e 


6. L, = JU*L 

1 W'> 


REMTECH INCORPORATED 


RTR 008-4 


7. * L = L^l + 0.447 Qy-1 )/|] M 

8. For y/L < 1 


Table 2 (Continued) 

5.22 + 4. 



41 (v% 


if = 1.- .189 exp (-5.3y/l.) -10.6 V exp[-y/L(9.87n» 

U d ' ' n=f 9.87n 2 +■ 28.1 


+ 28.1) /z ] 


9. For y/L < 1 

VS 

StRe. 



.40 


’A 



where for Pr * .7 

Nomencl ature 



X Q - B.L. running length to cavity 


Cj = 0.418 1 

W,H- See Fig. 1 


. C 2 = 0.435 

U e - Boundary layer (BL) edge velocity 


Cg = 1.87 : 

U d - Dividing streamline velocity 

10. 

T w -.Wall temperature 
P r - Prandtl number 

For y/L >1 

T aw - Adiabatic wall temperature 


v = U d 

Cp - B.L. edge heat capacity 
P e - B.L. edge density 

. H- 

For y/L > 1 

u e - B.L. edge viscosity 


x/ ,, . Va 

T e - B.L. edge temperature 


StRe L '*■ - Cj/fy/L -A) 

M e - B.L. edge Mach number 



U c - Viscosity evaluated at T c 


where for Pr * 0.7 

A - 0.178 

M d L 
Re. = ■ 

L u 

M c 

y - Distance along cavity reattachment 
wall from the external surface 

y - Ratio of specific heats 

12. 

| q - Local heating rate 

13. 

< ’ = s V C p (T d- T w ) 

j- - ’ - % ’ . . 


r The misprint which occurred in this equation In Ref. 15 was corrected by 
personal coiroiuni cations with the author. 
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Hg. ] - Schematic of Open Cavity Flow Structure 
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RAREFIED FLOW HEATING 


Aeroheatlng during hypersonic flight at rarefied conditions has become 
a problem of concern for Space Shuttle launch and entry. Mission require- 
ments dictate that during ascent and entry the External Tank (ET) will spend 
a significant amount of time In rarefied hypersonic flight. Consequently, 
several reports have been prepared, Refs. 3 to 7, which detail expected flow 
conditions, rarefied flight predictive methods, analysis of available data 
and computer program documentation for methods of computing rarefied flow 
thermal environments. The more salient features of this work are given here. 

The most significant contribution of Ref. 3 Is the correlation of rare- 
fied flow heating data In a manner which provides a rapid calculation proce- 
dure of rarefied flow heating for the design engineer. Flow regime switching 
criteria are given In Table 3. Relations for free molecular heating are 
given In Table 4, Rarefied flow heating relations and data are presented 
for spheres, cones, cylinders, flat plates and wedges In Tables 5 to 9 and 
Figs. 3 to 7. Wedge flow relations were developed subsequent to the work 
given in Ref. 3. In addition to the heating methods stated above. Ref. 3 
provides an alternate method of scaling wind tunnel rarefied stagnation point 
heating to flight conditions. 

The rarefied and free molecular flow heating prediction methods were 
programmed and Incorporated Into the MINIVER computer program. Ref. 6 pro- 
vides documentation of the changes made and Ref. 7 provides a program list- 
ing of the entire modified computer code. 

During the course of checking-out rarefied flow heating methods, It was 
found that a simple method was needed to compute the local recovery enthalpy. 
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Ref. 7 was completed in response to this need. A Newtonian and a concentric 
shock approach were shown to yield the same relationship for the recovery 
enthalpy. This relationship is given in Table 10. 

In order to attempt to broaden the rarefied flow heating data base, a 
review of Saturn V launch data was undertaken. The basic conclusions drawn 
from this study (Ref. 4) are: (1) The Saturn 501 heating rate gauge output 

in the rarefied flow regime is of the same magnitude as the zero bias error 
of the respective gauges, (2) Since some of the gauge output exceeds the 
maximum convective heating potential, gauge output cannot be attributed to 
convective heating and must be due to gauge inaccuracy and/or solar radia- 
tion input and (3) The uncertainties and inconsistencies in the data pre- 
cludes their use in establishing a basis for rarefied flow convective heat- 
ing. 


11 
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TABLE 3 

FLOW REGIME CRITERIA 


K, 

< 0.05 

Boundary layer flow 

^Re«> 



/ M* 

0.05 < 

< 3.0 

Rarefied flow 

yjKeco 




> 3.0 

Free molecular flow 
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. TABLE 4 


FREE MOLECULAR FLOW HEATING EQUATIONS* 


(1) 

N * 

P 

(particles/ unit volume) 

t 

U A 

co 

(2) 

S 

jZ* (molecular speed ratio) 

yp? 

(3) 

i 


NU„ 

ZtIvS ■. : 7“ ' 



(4) n “ S sin 0 

(5) n * 4* Je' 11 + fi T n(l + erf n)j 

v * ntanber of molecules striking a unit area per unit time 

where • 

' . A '• ' f 2 

erf n = — / e" A dx (error function) 

Rational approximation (o £ n < ®) from Ref. 5 
erf n ” 1. -(a^t + a^t 2 + a^vje" 1 * + e(n) 

t = -3- 

l+a 0 n 

a 0 * 0.47047 

a « 0.3480242 

a s -0.0958798 ' 

2 

a a 


0.7478556 
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* 

TABLE 4 (cont.) 

FREE MOLECULAR FLOW HEATING EQUATIONS 
(6) <f> = ~ e‘ n 


( 

/ . 


^ ■"» - 

fs J + jq n - «(, 

kT »| 

) 2 (Y-1 > W 

V y-V 



' * heating rate 
Nomenclature 

Uoo - Freestream velocity 

P - Freestream pressure 

T - Freestream temperature 

k - Boltzman's constant / 

Y - Freestream specific heat ratio 


♦From Qppenhelm (Ref. 17) 


H 


TABLE 5 
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f HEAT TRANSFER EQUATIONS ' 


Nomenclature 

T$ - Post normal shock temperature 
T^ - Wall temperature 
^ - Freestream flach number 
Y - Freestream specific heat ratio 
U f - Viscosity evaluated at T f and P g 
Hn - Total enthalpy 
Hy - Wall enthalpy 
Pa, " Freestream density 
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TABLE 6 

i • 

RAREFIED FLOW SHARP CONE 
HEAT TRANSFER EQUATIONS 


(1) T r * T w + (T 6 + T w )/2 - T 6 cos 2 0 c /3 



l°g iQ I a 1 (log 10 5i c ) 1 

• - ’ 1*0 


a 0 = -0.392510 : 

= -0.266308 
a* = -0.0598724 
(7) Heat Transfer 

R a p « u « c h v:v Nomenclature 

See Table 5 
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TABLE 7 

RAREFIED FLOW SHARP FLAT PLATE 
HEAT TRANSFER EQUATIONS* 



(5) C Hsi * (0.368VT o + 0.0684) [^(C*/ReJ Vi ] 


(6) 

CH 

Si 

= i-[l - tanh (0.91 log 10 g + 1.10)J 




for 3 < 0.1 


(7) 

Heat 

Transfer : . ; 



q s 

P «U„ C H (H. - H,) 

Nomercl ature 




See Table 5 


*From Shorensteln and Probsteln (Ref. 18) 
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TABLE 8 

RAREFIED FLOW YAWED CYLINDER 
STAGNATION LINE HP' TRANSFER EQUATIONS 


RTR 008 - 


(1) T r = i (T w + cos 1 


A) 


(2) Re 

(3) C* = 




Tr 


Re, 


A. 


U) l z * 

' \ Y«Moo 2 C*cosa 


Nomenclature 
See Table 5 


(5) S = fy/cosA 


(6) Heat Transfer Correlation 
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TABLE 9 

RAREFIED WEDGE FLOW HEAT TRANSFER 
■- .EQUATIONS 


(1) T r = T 6 /6 + 

(2) C* = ^ 

^aaT 


(3) = MJL 


WV 


R e«,x 


YdC C*cose 

(5) c - C H /COS0 


Nomenclature 

T^ - Post Normal Shock Temperature 
T w - Wall temperature 
T„ - Frees tream temperature 
X - Surface length 

- -Frees tream velocity 
Pw - Freestream viscosity 
P» '• Freestream density 

Y ro - Freestream specific heat ratio 
Moo - Freestream Mach no. 

Ho, - Freestream total enthalpy 

- Wall enthalpy 


(6) Heat Transfer correlation^ 0 > 10°) 

log 10 £ * E ^Oog,,^) 1 ; 

' '. 1=0 


a. =-0.370778 
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TA8LE 10 

A SIMPLE RECOVERY ENTHALPY RELATION 
FOR HYPERSONIC FLIGHT 


H r =* h w + (sin 2 6+ r cos 2 8) U w 2 /2gJ 


Nomenclature 


Hp.- Recovery enthalpy 

hoo - Freestream static enthalpy 

_i_ 

r - Recovery Factor = P r 2 (Laminar) 

i 



s P r * (Turbulent) 


P r - Prandtl number 


U* - Freestream velocity 
g - Gravitational constant 


0 - Mechanical equivalent of work 



1.000 


lifUiij 


ills; 


Free molecular 
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IBUtt 


( ■ i ■ .f ■ 
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Theoretical 

Cheng 
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Experimental Data 

Witt! iff & Wilson (1960) 
Ferri Zakkay & ting (1960) 
Hacker & Wilson (1963) 
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Fig. 3 Stagnation Point Heat Transfer Coefficient for Low Reynolds Number FT 
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Fig. 7 Low Density Wedge Flow Heat Transfer Data 
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Section 4 

IMPULSE MODEL RESEARCH 

The objective of this work was to develop a sound understanding of 
Impulse Base Flow Facility (IBFF) model and tan* wave dynamics from which 
future designs could be derived. To accomplish this objective, several 
work Items were completed. The performance of five IBFF model configura- 
tions were monitored and analyzed as they were used for other test programs 
(Ref. 8). Theoretical steady state H 2 -0 2 combustor properties were calcu- 
lated and presented In graphical form for design and/or analysis work by 
Ref. 11. A pretest plan (Refs. 9 and 10) was developed to obtain para- 
metric experimental data for design of impulse operated models and blast 
wave tank dynamics. The test was conducted and the data analyzed and pre- 
sented In Refs. 12 and 13. 

4.1 IBFF DATA MONITORING 

Facility model data was taken during routine testing operation of the 
IBFF. One hundred and twenty two (122) runs were made In which the combustor 
pressure. Initial combustor pressure rise rate, and type of combustion were 
recorded. This data was taken using two injectors and two nozzle configura- 
tions with which nozzle diaphragms of various thicknesses were used. 

From the data analyzed in Ref. 8 the following conclusions were drawn: 
Ignition Characteristics: 

• Ignition resulting in hot flows were experienced in all runs 
using the equivalent engine with a diaphragm and a pin hole. 

• It follows from the above conclusion that a glow plug is not 
necessary to obtain repeatable ignition of H 2 - 0 2 in some 
configurations. 
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• No correlation of pin hole size or diaphragm thickness with 
Ignition behavior was found. 

• Not using a nozzle diaphragm reduces ignition reliability 

• Insufficient data was available to quantitatively establish 
the difference between ignition characteristics of the large 
and small volume injectors in the equivalent engine. 

• The two engine geometry exhibited poor Ignition reliability. 
This configuration has a larger combustor volume due to the 
two nozzle entrance ports, 

t The two engine configuration is the only configuration studied 
which has produced an overpressure combustor pressure response. 

• The timing between H 2 and 0 2 arriving at the combustor may 
effect the percent hot flows obtained. 


Performance Characteristics 

• The chamber pressure for hot flows Is repeatable to within 
at most t 12% of the average chamber pressure for a given 
geometry, 

• Injector geometry and/or volume can significantly change 
the measured steady state chamber pressure. 

• The average chamber pressure for Initial overpressure runs 
Is statistically indistinguishable from the average chamber 
pressure for nominal hot flow runs. 

• The initial cold flow chamber pressure rise rate can be 
statistically distinguished from hot flow chamber pressure 
rise rates. 

• The hot flow initial chamber pressure rise rate is statis- 
tically the same for a 0.005 In. thick nozzle diaphragm as 
those runs without a diaphragm. 

• Runs which exhibit Initial overpressure characteristic also 
exhibit the largest scatter in the measured steady state 
chamber pressure (i.e. t 28.5% from the mean). 


4.2 CHAMBER PROPERTIES OF GASEOUS H 2 -0 2 COMBUSTION 

The purpose of this study {Ref. 11) was to produce curves showing the 
effects of Oxidizer - Fuel Ratio (0/F) and chamber pressure (P c ) on the 
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steady-state parameters of temperature, specific heat ratio (y), molecular 
weight and charactersitic velocity (C*) In the combustion chamber of an 
Impulse Base Flow Facility (IBFF) model. The curves presented in Ref. 11 
are Intended to provide the engineer with a rapid method for predicting 
combustion chamber conditions for gaseous fueled hydrogen-oxygen engines. 

4.3 IMPULSE MODEL TEST PROGRAM 

References 9 and 10 present the test plan developed to obtain parametric 
experimental design data for Impulse operated models and tank wave data. 
Specific test objectives which were accomplished are: 

1. Evaluation of model performance with different: 

• Combustor volumes (L* * 10,15,25,40 In) 

• Injector configurations (Triplet, doublet, concentric) 

i Charge tube pressures 

• Diaphragm effects (with or without) 

• Glow plug effects (with or without). 

2. Determine the behavior of the blast wave and plume Interference 
effects in the vacuum tank for variations In: 

• Tank pressure 

t Location In the tank 

Figure 8 presents a schematic of an impulse model installed in the tank 
bulkhead of the IBFF. The facility basically consists of two charge tubes 
(oxidizer and fuel) each containing a flow metering venturi, a diaphragm and 
cutter assembly, an Injector, a combustion chamber, a nozzle or nozzles and 
an altitude simulation chamber. Figure 9 presents an assembly drawing of one 
of the configurations tested. The triplet element injector used in these 
tests is shown as part of the assembly drawing. The other Injectors used 
are shown in Figs, 10 and 11. By replacing anular sections of the combustion 
chamber (see Fig. 9) the combustor L* was changed allowing four combustor 
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1 

volumes to be tested. The actual volumes and areas for the configurations 
tested are listed in Table 11. 

A detailed description of the test conducted and results obtained is 
given in Refs. 12 and 13. Reference 13 documents the test program with 
respect to the impulse model characteristics and Reference 12 documents the 
tank wave data. Reference 13 also gives a detailed discussion of model de- 
sign methods. Each major component of an impulse operated model is dis- 
cussed in terms of previous designs and recommended design methods. The 
I . • . : ... 

components considered are charge tube , venturi*, diaphragm cutters, diaphragms, 
injectors and combustors. 

In order to explain the data obtained through this test program and to 
provide a design tool for future work, an analytical time dependent combustor 
model was developed. Several simplifying assumptions were made to obtain a 
tractable model. The basic assumptions consisted of: 

1. One dimensional flow 

2. Constant combustor specific heat ratio 

3. Instantaneous heat release In the combustion chamber by the 
products present 

4. Temperature independent heat release 

5. Instantaneous perfect mixing within the combustion chamber 

6. Constant temperature throughout the manifolds 

7. Combustion chamber species were limited to H 2 , 0 2 and H 2 Q 

8. The nozzle throat diaphragm opens instantaneously at a specified 
pressure. 

Based on these assumptions, a time dependent math model for the transient 
behavior of the combustor was derived by writing mass, pressure, and energy 
balances. 

Balance equations were written for: 

1. H 2 mass balance in the volume preceding the injector face (H 2 mani- 
fold volume) 

2. 0 2 mass balance in the volume preceding the Injector face (0 2 mani- 
fold volume) 
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3. Pressure balance across H 2 and 0 2 Injector ports 

4. Combustion chamber mass balance 

5. Combustion chamber energy balance 

6. Combustion chamber species balance for H 2 , 0 2 and H 2 0. 

The mathematical form of these balances are given in Table 12 along with 
other equations required for complete definition of the combustor processes. 

A few comments are relevant to the equations given in Table 12. The 
injector pressure loss equations specify the velocity through the injector 

i 

ports leading to the combustion chamber. Choked flow through the ports was 
assumed until the chamber pressure rose to 52.8 percent of the respective 
manifold pressures. The energy balance equation was derived by specifying 
all input, output, and accumulation terms. The relation given by Eq. 8 was 
then obtained by subtracting Cp(T - 1080) times Eq. 7 from the original equa- 
tion. The effects of variable 0/F on Cp have been Included by including a 

multiple segment interpolation calculation of as a function of 0/F. 

■ ' 

In order to compare the results from the transient equation to steady 
state conditions, the steady state equations for the combustor were derived. 
The equations are given in Table 13 and consist of basically the same type 
as given in Table 12. 

The set of coupled time dependent equations given In Table 12 were 
solved numerically using a simple Euler's numerical integration method. 

This method was found to be satisfactory if sufficiently small step size 
was used. In all cases considered the transient solution obtained the steady 
state value if the integration was extended over a sufficient time Interval. 

A complete documentation of the computer code for the solution of the 
steady and transient equations is given in Ref. 13. 

The transient math model has been used to analyze the data obtained 
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during the test program. Fig. 12 presents a comparison of the theoretical 
transient response of the three main configurations studied experimentally. 

It can be noted that the doublet injector configuration obtained steady 
state more quickly than the other configurations. This was verified 
experimentally. Figure 13 shows a comparison of the theoretical method 
with data. The comparison is quite satisfactory for design work and is 
typical of the many cases studied in Ref. 13. 

Reference 12 presents the measurements made in the altitude simulation 
chamber or vacuum tank to determine the effects of the reflected blast wave 
on the nozzle exhaust plume for variations in tank pressure and tank loca- 
tion. The prime objective of the vacuum tank measurements was to determine 
the time interval for which steady state conditions existed in the exhaust 
flow before the reflected wave disturbed the flow. A secondary objective 
was to attempt photographic coverage of the plume during steady state condi- 
tions. To accomplish this objective color and Infrared film were used and the 
f-stop setting and development procedure were varied. Only marginal success 
was achieved. 

The reduced data presented in Ref. 12 for the primary test objective 
consists of different tables and plots of the time at which the Initial and 

• . i ■ 

reflected waves passed different probe positions. These plots were then used 
to construct a graphic rep resen tai ton of the wave patterns which existed at 
different time intervals after firing the model. The data was compiled for 
the three test altitudes, l.e. 240 Kft. , 280 Kft., and 320 Kft. All the 
wave data at a given altitude was averaged together and no distinction was 
made for different L*, Injector arrangements or diaphragm burst pressures. 
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The wave patterns given in Figs. 14, to 19 show the initial and reflected 
wave at different times for each of the three test altitudes, 240, 280, and 
320 Kft. For each altitude the Initial wave is shown at times 0.5, 1.0, 1.5, 
and 2.0 mill i sec. and the reflected wave is shown at times of 5.0, 7.5, 10. C. 
and 15.0 mini sec. For the Initial wave patterns constructed, the following 
items were noted. 

1. As the test altitude was increased the initial wave exhibited an 
expected broadening due to the expansion into a lower pressure. 

2. Likewise as the altitude Increased the wave speed Increased. At 
240 Kft. the wave was just reaching the tank end and side wall 

at 2.5 mlllisec. while at the 320 Kft the wave had almost reached 
the side wall at 1.5 mlllisec. 

3. The shape of the wave also changed with Increasing altitude. At 
the lower altitudes the wave front had a parabolic shape that 
changed to a wave front which was almost perpendicular to the tank 
centerline at 320 Kft. At each altitude the wave was found to be 
relatively flat and parallel with the nozzle exit for some distance 
from the exit. It appeared that this flatness was maintained for 

a longer time at the higher altitudes. 

4. The waves In the region which Is shown in the figures as extending 
from the lip of the nozzle toward the tank wall were found to move 
slowly In comparison to that movement in the frontal wave regions. 

In fact, there was strong evidence to support an argument that the 
phenomena observed was not that of a wave but simply the gas front 
expansion. This argument was supported by the fact that the phe- 
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nomena in most of this region was moving at a speed that was less 
than Mach 1 and that the thin film gauge responses did not exhibit 
the rather abrupt change associated with a strong blast wave. The 
movement in the region was so slow that the reflected wave from the 
tank side wall appeared to overtake the initial gas expansion. 

The reflected wave patterns are of even greater significance to short 
duration rocket motor tests at simulated altitudes because the arrival of the 

I 

reflected wave within the plume or model base region destroys the altitude 
simulation and signifies the end of the test time. It was unfortunate then 
that the reflected wave patterns were more difficult to establish and identify 
than the initial wave. This was because at the time of the passage of the 
reflected wave, the thin film gauge was still responding to the passage of 
the initial wave and the gas following it. Although this In Itself made 
the identification of the passage of the second wave difficult, in addition, 
the reflected wave was weaker than the Initial wave and therefore did not 
produce a pronounced change In the gauge response. For these reasons the 
reflected wave data was found to be more scattered than the initial wave 
data. For the reflected wave patterns constructed, the following items were 
noted. 

1. The reflected wave and/or gas movement region which was likely to 
first cause steady flow breakdown was found to occur from the tank 
side walls and to form a "U" shaped pattern which generally moved 
along the plume boundary. 

2. Although the reflected wave movement pattern differed somewhat 
with altitude, it was found to arrive within the nozzle base region 
at approximately the same time. Average arrival times at the tank 
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bucket corner were 11.8, 11.0, and 9.8 msec for 240, 280, and 
320 Kft. altitudes respectively. 

3. As noted earlier in the discussion of the initial wave, the region 
from the nozzle lip to the tank side wall has a very weak initial 
wave and is probably only the gas expansion in that region. For 
this reason it is likely that no well defined reflected wave comes 

directly off the side wall. The sketched reflected wave patterns 

I 

are most likely a combination of wave phenomena and gas front ex- 
pansion. The ”U U shaped wave, even though it is produced within 
a region where the reflected wave is very weak, moves rapidly back 
to the nozzle region because it has very little or no gas velocity 
head to move against, thereby achieving a rather quick relative 
movement. 

4. There was some evidence of the existence of a third wave phenomena. 
However, the data obtained were Inconclusive, so no definition of 
its wave pattern was made. This wave was suspected to come from the 
region where the tank far end wall joins the tank side wall. The 
wave movement would be along the tank side wall and behind the wave 
reflected from the tank side wall. 
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12 Element 
Triplet 


H * 

Venturi Diffuser + Diaphragm 
Bulqe Volume (In. 3 ) 

0.951 

0.914 

‘Injector Housing - Diaphragm 
Bulqe Volume (in. 3 ) 

0.474 

1.150 

Injector Volume (In. 3 ) 

1.958 

4.493 


2.432 

mm 

Injector Face Element 
Area (In. 2 ) 

0.782 

0.372 

Combustion Volume 
In Injector (in. 3 ) 

2.946 



. to 
Throat 

Nozzle and Adapter 

3.54 

6.12 

Nominal L* = 10 Comb. Section 

4.47 

4.47 

Total for L* = 10 

10.96 

13.54 

Actual L* (In.) for L* - 10# 

12.4 

Nominal 1* « 15 Comb. Section 

9.84 

9.84 

Total for L* = 15 

16.33 

18.91 

Actual L* (In.) for L* = 15 

17.3 

Nominal L* * 25 Comb. Section 

20.83 

20.83 

Total for L* = 25 

27.32 

29.90 

Actual L* (in.) for L* = 25 

27.4 

Nominal L* - 40 Comb. Section 

37.20 

. 37.20 

Total for L* = 40 

43.69 

46.27 

Actual L* In.) for L* = 40 

42.5 


# A* = 1.091 In. 2 
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6 Element 
Triplet 

Concentric 

Doublet 






H 

2 

0.951 

0.914 

0.951 

0.914 

0.951 

0.914 

0.474 

1.150 

0.474 

1.150 

0.474 

1.150 

1.542 

4.493 

IMBfrll 

8.253 

0.621 

3.054 

2.016 

5.643 

3.077 

9.403 

1.095 

4.204 

0.391 

0.186 

0.964 

1.374 

0.230 

0.230 

2.946 

2.936 

6.193 

Vol 

Diaph. 

to 

Throat 


to 

Throat 

Vol 

Diaph. 

to 

Throat 

3.54 

6.12 

3.54 

6.12 

3.54 

6.12 

4.47 

4.47 

4.47 

4.47 

4.47 

4.47 

10.96 

13.54 

10.95 

13.53 

14.20 

16.78 

12.4 

12.4 

15.4 

9.84 

9.84 

9.84 

9.84 

9.84 

9.84 

16.33 

18.91 

16.32 

18.90 

19.57 

22.15 

1 17.3 

17.3 

20.3 

20.83 

20.83 

20.83 

20.83 

20.83 

20.83 

27.32 

29.90 

27.31 

29.89 

30.56 

33.14 

1 27.4 

2 

7.4 

30.4 

37.20 

37.20 

37.20 

37.20 

37.20 

37.20 

43.69 

46.27 

43.68 

46.26 

46.93 

49.51 

4 

2.5 

42.4 

45.4 
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TRANSIENT EQUATIONS FOR GASEOUS PROPELLANT 
IMPULSE OPERATED COMBUSTORS 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

(6) 
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TABLE 12 (Cont. 1} 


Combustor Mass Balance 

dp_c 
dt 

where 


Vc fje * (pUA) 0 + (pUA) H - p c KC c A*D 


(7) 


0*0 Before the nozzle diaphragm burst 

/ ’ ■ ’ • 

* 1 After the nozzle diaphragm burst 

(Y+l }/2(Y-l ) 

Y 3 1.135 assumed for O/F = 6.0 




Combustor Energy Balance 
dT 


p c V c C p ^ = f Hr w - [(pUA) h + (pUA)J Cp(T-1080) 

-126(pUA) 0 - 1903(pUA) H + j 

where H = heat of reaction (5730 Btu/lbm of H 2 0) 
Combustor Species Balance 


0 2 - Balance 
dW. 


W “ (p ua > 0 - 5f p c A ‘ C c® - r o 


where 


r 0 = <P UA ) 0 1f W 0 /W H 5 8 

- |% 1f V W H > 8 


( 8 ) 


(9) 


Ho - Balance 


dW H 3 (pUA) H - % 


dt 


— p r A*C c KD - r H 

Wy c C M 


00 } 
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where 

r H = ( pUA ) H if W Q /W H > 8 


l 

9 


dW 

Z* If W 0 /W H £ 8 
dt 0 H 


H 2 0 - Balance 


dW, , / w, , 

— = r w " , PrA*C c KD 

dt w W T c c 


where 


r w = | dW„ 1f W(>/Wh ^ 8 

dt 


* 9 If W 0 /W„ > 8 
dt ° H 



Combustor Mass 

«r " w h + w o + W w 
Combustor OF 

(° F ) C = (W 0 + 8W w /9)/(W h + W J9) 


Combustor Molecular Weight 
M c = (W H + Wj/OyZ + W^/18) If W 0 /W K < 8 

M C = (W o + w w )/(y 32 + v 18 ) if v w h > 8 
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02 ) 

03 ) 


04 ) 


05 ) 
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STEADY STATE EQUATIONS FOR GASEOUS PROPELLANT 
IMPULSE OPERATED COMBUSTORS 


Combustor OF 

(0F) C = m 0 /ra H (1 ) 

Combustor Molecular Weight 

M c •- 2(0F) C + 2 If (0F) C 5 8 (2) 

M c - 32 ( (0F) C + 1 }/ ( (OF ) c + 8) If (OF) c > 8 (3) 

Combustor Species 
Before combustion (mass ratios) 

W 0 /Wr * (0F) c /((0F) c + 1} (5) 

V W T ’ 1 - V“ T ' 1 (6) 

Before combustion (mole to mass ratios) 

Nh/Wt - \ V W T (?) 

V W T = -l2 W o /W T (8) 

Post combustion (mole to mass ratio of H 2 0) 

yW T - 2 N o /M t If (0F) c ^ 8 (9) 

> N q /W t If (0F) C >8 (10) 
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TABLE 13 (Cont* 1) 

Post combustion {mass ratio) 

yUj = 18 yW T (U) 

Combustor Energy 

T «' H(VW T )/C p (12) 

Speed of Sound 

C c s ( y c RT9 c /M c) 03) 

Combustor Mass Balance 

P c = (■„ + m H )/(KC c A*) (14) 

wnere / , \ (Y+1)/2(Y-1) 

‘ ' (- ) 

Combustor Pressure 


P c » PcRT/M c 05) 

Injector Face Densities 

Poc = W RT o 06) 

. . ’ ’ 1 . • ' 

P H c = P c h h /rt h 07) 

Injector Velocities 

u o = V(Voc ) 08) 

U H • Whc* 09) 

■ 40 ’ ' 
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TABLE 13 (Cont. 2) 



P o “ P c + P M U o p o/9 c 
P H = P c + PHc U H F H /9 c 


(20) 
( 21 ) 




Nozzle 
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v di i au j«_ Lwi f y ui 

Combustor Wall 



Pressure Gauge Port 


H 2 Upstream and Dcv/nstream 

Pressure Gauge Locations (out of plane) 


Fig. 9 Combustor Assembly, L* = 10 In., Triplet Injector 
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Time (msec) 

Fig. 12 Theoretical Start-up Processes of 8% 0-2 Combustor with the Doublet and Triplet and Concentric 
Injectors for L* = 25 and P CB^^ P C ^steady = °* 956 


REMTECH INCORPORATED RJR 008-4 


Combustor Pressure (psla) 


CB = 600 psla 
= 550 psla 
P CB * 500 psla 


iSasBiisHi 


1 1st 


Doublet Injector 
L* = 25 in. 

3(0.005) In. Nozzle Diaphragm 


• - 58/0 
A- 59/1 
■ - 60/0 




114 I I t H MU I 1-14 1 t 


Time (msec) 

Fig, 13 Comparison of Experimental and Theoretical Combustor Pressure Histories for the Doublet Injector 
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Section 5 
CONCLUSIONS 

The work conducted under this contract included several technical areas. 
These areas have been summarized In the preceding section and fall within 
the following groupings. 

1. Cavity Heating (Ref. 1 and 2) 

, 2. Rarefied Flow Heating (Ref. 3 to 7) 

3. Impulse Model Research (8 to 13) 

Specific documentation has been given In each of these groupings as indicated 
by the noted references. In each of the three groupings either methodology 
Improvements were achieved or experimental data has been obtained which 
provides a more complete understanding of the areas studied. 
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